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§ ! Abstract 

.' In strongly interacting electron systems with low density and at low temperature the thermo- 

r^ ■ dynamic density of states is negative. It creates difficulties with understanding of the Einstein 

f^^ . relation between conductivity and diffusion coefficient. Using the expression for electrochemical 

potential that takes into account the long range part of the Coulomb interaction it is shown that at 

IJ ' negative density of states Einstein relation gives a negative sign of the diffusion coefficient D, but 

U 

'jrt '• unipolar relaxation of inhomogeneous electron density is not described by the diffusion equation. 

The relaxation goes much faster due to electric forces caused by electron density and by neutral- 
P^ ■ izing background. Diffusion coefficient is irrelevant in this case and it is not necessarily positive 

o: 

O . because process of diffusion does not contribute to the positive production of entropy. In the case 

«vj , of bipolar diffusion negative D results in a global absolute instability that leads to formation of 

(j\^ ' neutral excitons. Graphene is considered as an example of a system, where the density relaxation 

Q> ■ is expected to be due to electric force rather than diffusion. It may also have a negative density of 

o 

OO 

o 



under this condition there is no thermodynamic limitation on the sign of D. It happens because the 



states. 



PACS numbers: 71.27. +a,73.50.-h 
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I. INTRODUCTION 

The idea of the Einstein relation was put forward by Einstein^ and Smoluchowski^ in 
1905-1906. Both scientists considered the Brownian motion in the presence of gravitational 
force. The result is the relation between mobility u in the field and diffusion coefficient D. 
In case of electric field and particles with the charge e it has a form 

eD = Tu, (1) 

where T is the temperature in energy units. The main idea was equivalence of an exter- 
nal force and the density gradient. Of course, both Einstein and Smoluchowski did not 
care about negligible mutual gravitational or any other small interactions of the Brownian 
particles. 

The formulation of the Einstein relation for electrons is based upon electrochemical po- 
tential, the thermodynamic function that, like temperature and pressure, should be the same 
at all points of the system in the equilibrium state. The usual arguments are as follows. If an 
external potential ip is applied to the system, the condition of thermodynamic equilibrium 
reads 

Eec = /^(^) + eip = Const, (2) 

where //(n) is the chemical potential as a function of inhomogeneous electron density n. 
In the equilibrium both n and ip are function of coordinates while E^c is constant. The 
temperature T should also be constant. Therefore, the electrical current density j at constant 
T can be written in a form- 

i = --VEec = aE-DVen, (3) 

where a is conductivity and E = —Vip. Then one gets relation connecting a and D 

which is also called Einstein relation. For the Boltzman gas d^/dn = T/n and one gets 
Eq-dID if 0" = enu. It looks like derivation of Eq. (jl]) is independent of the properties of the 
system and this equation can be consider as general thermodynamic law. 

A simple observation shows however that in the case of non-ideal electron gas the Einstein 
relation needs some comments. We discuss an electron gas on the positive background at 



low temperatures and low densities when dimensionless parameter r^ is not very small. Here 
r^ = 3/{ATTna^) for 3-d case and r^ = l/7m2a^, where n and n2 are 3- and 2-dimensional 
electron densities respectively and a^ = ffK/me^ is the Bohr radius, m is an effective 
electronic mass, k, is an effective permittivity. 

The problems of dynamic screening and diffusion in slightly non-ideal electron gas (r^ << 
1) with electron-electron interaction were considered in details about 20 years ago (See 
Ref. [^1^1^]) In this case the thermodynamic density of states is large and positive. I concentrate 
here on the strongly non-ideal case r^ > 1. 

An electron gas on the positive background at low temperatures and low densities has 
energy E of the order of —e^n^^'^'N/K, where c? = 2, 3 is the space dimensionality and n is the 
density per area or volume respectively, A^ is total number of electrons. Then yU ~ —e^n^^'^/K 
and E, /i, and dfi/dn are negative^i^. The first experimental confirmation of this idea was 
done by Kravchenko et a^^, but direct quantitative study of this effect was performed by 
Eisenstein et a^^. 

The derivative dn/dn is proportional to the reciprocal compressibility of the electron 
gas. Note that compressibility has to be positive due to the thermodynamical condition of 
stability. However, this principle cannot be applied to the charged systems, like electron gas, 
because part of their energy is outside the system in a form of the energy of electric field. 
On the other hand, in the case of a neutral electron-hole plasma, the situation of negative 
compressibility can arise leading to collapse of the system. Such a situation is considered at 
the end of Sec. HI. 

It follows from Eq. Hlthat if dfi/dn is negative, diffusion coefficient D and conductivity a 
have opposite signs. This observation needs an explanation because near the thermodynamic 
equilibrium both of them have to be positive to provide positive entropy production due to 
the Joule heat and due to the relaxation of inhomogeneous density. 

II. ELECTROCHEMICAL POTENTIAL AND STATIC SCREENING 

To resolve this contradiction one should include the long-range part of the Coulomb 
potential created by inhomogeneous electron gas into the function E^c in Eq. ([2]). This 
contribution is a functional of n{r). 

To find Eec taking into account electron-electron interaction one should minimize the 



Helinholtz energy F with respect to electron density n(r) at a given value of T and A^. For 
low T one gets 

e^ f f n'(r)n'(r')d^rd^r' f ^ , ,, . 



/e„V)«'r-^„/„'(r)d3, (5) 



where / is the Helmholtz energy density of a homogeneous electron system that results from 
the interaction in a neutral system, like the Wigner crystal or "Wigner liquid". Since this 
interaction comes mainly from the nearest neighbors and n(r) is a smooth function, one 
may assume that both / and chemical potential /i = df /dn are local functions of n(r) We 
assume also that n(r) = n + n'(r), where n is average density and n' <^ n. 
Minimization of this expression with respect to n' gives the equation 

„ , , , du , e'^ f n'(r')d^r' . , 

dn K J |r — r I 

It differs from Eq. ([2]) by the potential of electrons in the right hand side. Note that 
this potential is due to the violation of neutrality in a scale much larger than the average 
distance between electrons. To check this equation we consider thermodynamic equilibrium 
and find equations for the Thomas-Fermi static screening in 3- and 2-dimensional cases. 
Since E^c is independent of r in thermodynamic equilibrium one may take E^c — IJ'{n) as a 
reference point for the total potential ip defined as 

(7) 

n J I r — 1"' I 
It follows from Eq. ([6]) that 
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The Poisson equation has a form 

where pext is density of external charge. Using Eq. ([H]) one gets final equation for the 3-d 
linear screening 

V72 2 47rpexi tir\\ 

V v? = -q^^ . (10) 

Here 

9 47re^ dn , ^ 
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is the reciprocal 3-diniensional screening radius. 

Consider now a thin layer (x-y plane) with 2d electron gas separating two media with 
dielectric constants ki and /t2- In this case one should substitute n ^ n2S{z) and k ^ R = 
{ki + fi;2)/2. The results is^ 

W = -q2^5{z) - ^, (12) 

where 

2'Ke^dn2 ,„. 

92 = ^— -1-. (13) 

K djjL 

It is important that Eqs. fITOl) . flT2|) are applicable only if the screening is linear {n' <C 
n)r^. There is another serious problem of applicability the Thomas-Fermi approximation 
in the case of the negative density of states. Indeed, the dielectric permittivity in this 
approximation has a form 

<<i) = <^ - §) (14) 



g2 



in 3-d case and 



e(g) = «:(! - M) (15) 

in 2-d case. In both cases it has roots at q = \q3\,\q2\- The expression for the screened 

potential (p has a form 

/ N /' y?o(q)exp(iq-r)c?q 

^''' = / ^) • ''" 

where (po is a bare potential. Thus, the roots of e transform into the first order poles without 
any reasonable way of the detour. Such a detour follows from the casuality for the cj-plane 
but not for the q-plane. Moreover, the electrostatic potential should be real and one cannot 
add a small imaginary part in the denominator. Therefore I think that the poles do not 
have any physical sense. 

The reason is that negative sign of the density of states appears when q^, q2 are of the order 
of average distance between electrons f. At such distances the very concept of macroscopic 
field does not have sense. However, if the bare potential has only harmonics with q <C 
I gal, \q2\, the Eqs. (11011121) have a sense. Consider, for example, the screening of the positive 
charge Z at a distance Zq from the plane with 2-d gas(plane z = 0. The solution of Eq. (fT2|) 
has a {otidf^ 



where p is a polar radius in the plane z = 0. Suppose that \q2\z0 ^ 1- Now the contri- 
bution to integral Eq. flT7|) from g ~ \q2\ is exponentially small and one can ignore q in the 
denominator. Then 

"(") = , f:° .,3/. - (18) 

Kq2 [z^ + p^) ' 
Note that at ^2 < a positive charge creates a small negative potential in the plane with 

electrons. That is what I call " overscreening" . 

Extra electron density, as calculated from Eq. ([8]) is 

, ^ Zq . . 

It is negative and independent of the sign of g2- One can see that the total charge 

en2Trpdp = -Z (20) 

'0 

Due to geometry of the problem electric field is zero below the plane with electrons. As 

follows from Eq. ([8]) , the signs of charge density and potential are opposite if the density of 

states is negative. 

For the case of two such planes (double quantum well structure) Luryi-^^ has predicted 
a small penetration of electric field through the first plane. He has considered the case of 
positive density of states. Then the small penetrating field between two planes has the same 
direction as the incident field. 

Eisenstein at ai— studied this effect experimentally and found out that at negative density 
of states the propagating field is opposite to the incident field and this is also a result of the 
overscreening (see the quantitative theory in Ref.— >^>^). 

Negative density of states was also used^ for the explanation of magnetocapacitance data 
by Smith at alr^. 

III. CONDUCTIVITY VERSUS DIFFUSION 

Now I come back to the problem of the negative diffusion. If the system is not in 
equilibrium the electric current can be written in the same form as Eq. ([3]) 

j = --WE,,. (21) 

e 



Using Eq. ([6]) one gets 



j = aE - DVen' - a-V [ V^J^^}^_ (22) 

•^ K / r-r' ^ ' 



Here D is connected to a by the Einstein relation Eq. (j4]). Considering relaxation of the 
charge density one can ignore external field E. The relaxation is described by the continuity 
equation 

^-V, (23) 

or 

dien) f I dn^^, ,, A'Ken'\ , ,, 

The ratio R of the first (diffusion) term in the right hand side to the second (field) term is 
R = (ggL^)^^, where L^^ = V^n' /n' is the characteristic size of the extra charge and q^ is 
given by Eq. ( ITTI) . If electron gas is non-ideal, qs ~ 1/f, where f is the average distance 
between electrons. However, the very concept of diffusion equation is valid at L ^ f. This 
means that for the non-ideal gas |i?| <^ 1 and the diffusion term in Eq. ([21]) should be 
ignored. Then the equation has a simple solution 

n'(r, t) = n'(r, 0) exp -(Vtm), (25) 

where tm = /«/(47rcr) is well-known Maxwell's time. Coefficient D does not enter in this 
case in the entropy production and it does not have a physical sense. Thus in 3-dimensional 
non-ideal electron gas negative dfi/dn does not create any contradiction with the Einstein 
relation. 

In the 3d gas of high density /i ~ n'^^^ and R ~ (f/LY/rg with r^ < 1. In this case R 
might be large and diffusion is possible. However dfi/dn > 0, and D > 0. 

Now we consider the relaxation of the charge density in 2-dimensional case. Instead of 
Eq. ([24]) one gets 

d{en2) A dfi 2f /n 

_e_^,fn'{r')d-r'_ (26) 

K J |r — r'l 

Here n2, cr2 and V are 2-dimensional density, conductivity, and 2-dimensional gradient re- 
spectively. To consider the ratio R2 of the first (diffusion) term to the second (field) term it 



is convenient to make the Fourier transformation. Then one gets 

where ng is the Fourier transformation of n'g. 

Now we find that the ratio of the first ( diffusion) term in the right hand side of Eq. 
fl27|) to the second (field) term R2 = q/q2, where g2 is given by Eq. flT3|) . Similar to the 
3d case in the non-ideal gas \q2\ ~ 1/f and diffusion should be ignored. Then we get the 
Dyakonov-Furman equation^ 

^ = -vqug, (28) 

where velocity v = 27rcr2//€. The physical meaning of this equation is that extra density of 
electrons localized initially at some spot propagates in all directions with velocity v conserv- 
ing the total amount of extra electrons. Of course, this way of relaxation is more efficient 
than diffusion (random walk), because r ~ f t while r ~ vDt in the case of diffusion. Thus, 
diffusion coefficient D is irrelevant and negative dfi/dn does not create any contradiction 
with the Einstein relation In a high density electron gas R2 = qf/rg and diffusion mechanism 
is possible. In this case dfi/dn > and D > 0. 

One can consider this problem from a different point of view. In both 3d and 2d cases the 
negative diffusion coefficient D appears in the term with the highest derivative that leads 
to the absolute instability even if D is smalP^. Consider, for example Eq. (!24|) for 3d case. 
After the Fourier transformation the solution for the charge density p = en' can be written 
in a form 

Pg = pI exp (-^ - DqHj , (29) 

where D is given by the Einstein relation Eq. (jlj). One can see that at D < solution 
increases with time exponentially for harmonics with qf > 1. 

The physical explanation is as follows. The Eqs. (l24|26p contain average distance between 
electrons f. So they contain information that the charged liquid has a discreet electronic 
structure. This information comes from the negative density of states which originates 
from the interaction of the separate electrons. That is why macroscopic equations become 
unstable at small spacial harmonics. The message is that n{r) is rather a set of 5-functions 
than a continuous function. The instability is absent if D is positive. 

The instability of small spatial harmonics at small negative D does not affect larger 
harmonics because Eqs. (l24|26p are linear. Due to the linearity different harmonics are in- 
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dependent and transformation of energy from small spacial harmonics to large harmonics is 
forbidden (cp. phenomenon of turbulence in non-linear hydrodynamics where the transfor- 
mation of energy is not forbidden, but the instability is initiated by large harmonics). 

Therefore, I think that at small D approximation D = Q that gives Eqs. fl25ll28p is correct. 

One should note that the problem of the non-physical roots of electric permittivity dis- 
cussed in the previous section is of the same nature. 

Before we discussed the unipolar diffusion. Consider the simplest case of the ambipolar 
diffusion assuming that at t = the densities of electrons and holes are equal in some finite 
region of space and are zero otherwise. Moreover we assume that the local macroscopic 
charge density p(r, t) = and a recombination of carriers is very slow. In this case Eq. ^ 
describes the electron- hole system in quasi-equilibrium. At large r^ one gets E, fi, dfi/dn < 
but the last term in Eq. ([6]) is absent. So the smearing of the density of particles is described 
by the equation of diffusion at all r^, but at small density (r^ > 1) coefficient D < 0. Then 
the absolute instability takes place for all harmonics that means a collapse of the system. 
Thus the electron-hole "Wigner liquid" and crystal are unstable. 

This result is very transparent. It happens because negative /i just means that the energy 
of the system decreases with increasing density. In bipolar case neutrality is provided by 
the particles and we do not consider any background. Thus the instability is a result of the 
negative compressibility in a neutral system. At large enough r^ these particles are classical, 
and the absence of the mechanical equilibrium follows also from the Earnshaw theorem. In 
reality quantum mechanics becomes more important with increasing density. As a result 
the excitons are formed. These neutral particles have a positive diffusion coefficient Da and 
their density smears with time through all available space. This process is described by a 
regular diffusion equation. In the case of optical excitation the carriers may appear in the 
form of the excitons from the very beginning 

For the coefficient of the ambipolar diffusion Da a textbook equation^ 

Da = —^ (30) 

"" D, + Dh ^ ^ 

is often used, where D^^h are diffusion coefficients of electron and holes in unipolar case. 
As follows from the previous discussion, one should be careful with this equation because 
for the non-ideal electron (or hole) gas these unipolar coefficients might be negative and 
meaningless. It happens because in unipolar case there is a deviation from neutrality that 



creates electric field, while in bipolar case the system is neutral. In this case Eq. (130|) does 
not work and one should calculate Da in a different way as a diffusion of the exciton. 

In the recent paper by Zhao^^^ the experimental results for the ambipolar diffusion in 
silicon-on-insulator system are compared with Eq. fl5U]) . At high temperatures a good 
agreement is found while at low temperatures the observed values of Da are 6-7 times less. 
The previously reported values^i show similar temperature dependence. 

The author's explanation is that coefficients D^.h are taken for the bulk silicon using 
Einstein relation and they might be larger than in the film at low temperatures. However, 
the reason discussed above cannot be excluded. 

IV. GRAPHENE AS A POSSIBLE EXAMPLE OF A NON-IDEAL ELECTRON 

SYSTEM 

It is interesting to discuss the single layer graphene as an example of the system with 
non-ideal electron gas. Graphene is a gapless material with the linear spectrum of electrons 
and holes near the Dirac point. Due to some reasons, that are not quite clear now, the 
velocity v of electrons and holes in equation e = ±pf is of the order of e^ /h. It follows 
that at any Fermi energy inside this linear spectrum electron gas in graphene is non-ideal 
in a sense mentioned above: the absolute value of the chemical potential is of the order of 
interaction energy e^n^^"^. It means that unipolar density relaxation in this system should 
be described by the Dyakonov-Furman equation rather than by diffusion equation. 

However, without magnetic field the electron gas in graphene is marginally non-ideal. It 
cannot be classical, like an electron gas of a low density with quadratic spectrum. The 
marginal situation makes theoretical calculations very difficult. Nevertheless, it is accepted 
that the Wigner crystal in single layer graphene is absent without magnetic field2^»2^. The 
sign of dfi/dn is also an interesting question but very difficult for theoretical study. Recently 
tunneling microscopy experiment has been done by Martin et al.—. They claim that their 
measurement give the thermodynamic density of states and that it is positive. The last 
statement might be a result of disorder. 
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V. CONCLUSION 

Finally I argue that the negative sign of diffusion coefficient that follows from the Einstein 
relation at negative density of states does not lead to any contradiction because diffusion 
coefficient is irrelevant for the unipolar transport under this condition. The sign of the 
diffusion coefficient in this case should not be definitely positive because the diffusion is not 
the main source of the entropy production. In bipolar situation negative diffusion means 
the collapse of the system and formation of neutral excitons. 

I am grateful to Boris Shklovskii and Yoseph Imry for important discussion. I am espe- 
cially indebted to David Khmelnitskii and Emmanuel Rashba for multiple discussions and 
criticism. 



Electronic address: |efros@physics.utah.edu| 



^ A. Einstein, Annalen der Physik 17, 549 (1905). 
^ M. von Smoluchowsky, Annalen der Physik 21, 756 (1906). 

^ L. D. Landau and E. Lifshitz, Electrodynamics of Continuous Media (Butterworth-Heinenann, 
1984), chapter III. 

4 B. L. Altshuler, A. G. Aronov, and P. A. Lee, Phys. Rev. Lett 44, 1288 (1980). 

5 A. Y. Zyuzin, JETP Lett. 33, 360 (1981). 

^ B. L. Altshuler and A. G. Aronov, Electron- Electron Interaction in Disordered Systems (ed. by 
A. L. Efros and M. Pollak) (North-Holland, Amsterdam, 1985), p. 37. 

'^ B. I. Shklovskii and A. L. Efros, JETP Lett. 44, 669 (1987). 

^ M. S. Bello, E. I. Levin, B. I. Shklovskii, and A. L. Efros, Sov. Phys JETP 53, 822 (1981). 

9 S. V. Kravchenko, V. M. Pudalov, and S. G. Semenchinsky Phys. Lett. A 141, 71 (1989). 
^'^ S. V. Kravchenko, D. A. Rinberg, S. G. Semenchinsky, and V. M. Pudalov, Phys. Rev. B 42, 
3741 (1990). 

11 J. P. Eisenstein, L. N. Pfeiffer, and K. West, Phys. Rev. Lett 68, 674 (1992). 

12 J. p. Eisenstein, L. N. Pfeiffer, and K. West, Phys. Rev. B 50, 1760 (1994). 

13 T. Ando, A. B.Fowler, and F. Stern, Rev. Mod. Phys. 54, 437 (1982). 
1^ A. L. Efros, Phys. Rev B 45, 11354 (1992). 



11 



15 S. Luryi, Appl. Phys. Lett. 52, 501 (1988). 

16 A. L. Efros, F. G. Pikus, and V. G. Burnett, Solid State Comm. 84, 91 (1992). 
i'^ F. G. Pikus and A. L. Efros, Phys. Rev. B 47, 16395 (1993). 

1^ A. L. Efros, Phys. Rev. B 45, 11354 (1992). 

19 T. P. Smith, W. I. Wang, and P. J. Stiles, Phys. Rev. B 34, 2995 (1986). 

20 M. I. Dyakonov and A. S. Furman, Sov. Phys. JETP 65, 574 (1987). 

21 I am grateful to E. I. Rashba for this comment. 

22 K. Seeger, Semiconductor Physics. An Introduction (Springer, 1999), p. 124. 

23 H. Zhao, Appl. Phys. Lett. 92, 112104 (2008). 

24 M. Roshng, H. Bleichner, P. Jonsson, and E. Nordlander, Appl. Phys. Lett. 76, 2855 (1994). 

25 H. P. Dahal, T. O. Wehling, K. S. Bedell, J.-X. Zhu, and A. V.Balatsky, arXivxond- 
mat/0706.1689. 

26 R. Cote, J.-F. Jobidon, and H. A. Fertig, arXiv:cond-mat/0806.0573. 

2^ J. Martin, N. Akrman, G. Ulbricht, T. Lohmann, J. H. Smet, K. von Klitzing, and A. Yacoby, 
Nature Physics 4, 144 (2008). 



12 



